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Abstract. We study the topology of quasiperiodic solutions of the vortex filament equation 
in a neighborhood of multiply covered circles. We construct these solutions by means of a 
sequence of isoperiodic deformations, at each step of which a real double point is "unpinched" 
to produce a new pair of branch points and therefore a solution of higher genus. We prove 
that every step in this process corresponds to a cabling operation on the previous curve, 
and we provide a labelling scheme that matches the deformation data with the knot type of 
the resulting filament. 



1. Introduction 

In this sequel to [6] , we continue our study of the role of integrability for periodic solutions 
of the Vortex Filament Equation (also known as Localized Induction Equation) 

dt dx dx 2 ' 

a model of the self-induced dynamics of a vortex line in a Eulerian fluid, described in terms 
of the evolution of the position vector ~y(x, t) of a space curve parametrized by arclength x. 
Hasimoto's transformation [16] 



q = | k exp y J T dsj , (2) 

given in terms of the curvature k and torsion r of 7, maps the Vortex Filament Equation 
(VFE) to the focussing cubic nonlinear Schrodinger (NLS) equation 

iq t + q xx + 2\q\ 2 q = 0, (3) 

with the NLS potential q defined up to multiplication by an arbitrary constant phase fac- 
tor. The integrability of the NLS equation [T2], [32] implies that the VFE inherits many of 
the properties of the integrable equation, including a family of global geometric invariants 
(conserved quantities) [22J, a bihamiltonian structure [31 [221 123 [2H] , an d special solutions: 
solitons, finite-gap solutions, and their homoclinic orbits [U El [SI ED] ■ 

Periodic boundary conditions for the VFE give rise to closed curves; of those, the class 
of vortex filaments corresponding to periodic and quasi-periodic finite-gap NLS potentials 
provides contains striking examples of curves exhibiting special geometric features (such as 
symmetry and periodic planarity) and interesting topology. Our previous article [6] focused 
on characterizing geometric properties of finite-gap VFE solutions such as closure, symme- 
tries, self-intersection, and planarity in terms of the Floquet spectrum of associated finite-gap 
NLS potentials. In contrast, the current work concerns the topological information contained 
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in the algebro-geometric data of a closed VFE solution associated with a periodic finite-gap 
NLS potential. 

Before describing our approach to this problem, we will briefly introduce some of the main 
tools and results used in the paper. 

The NLS linear system and the curve reconstruction formula. The AKNS linear 
system for (J3]) consists of a pair of first-order linear systems |12j : an eigenvalue problem 

Ci<t> = A<£, (4) 

and an evolution equation 

4> t = £i4> (5) 

for the complex vector- valued eigenfunction <fi. The solvability or "zero curvature" condition 
of the AKNS system is the NLS equation Q. Expressed in terms of the Pauli matrix 
(\ 

The coefficients of the linear operators C\ and £2 depend on x and t through the complex- 
valued function q, and on the spectral parameter A G C. 

The inverse of the Hasimoto map (i.e., the reconstruction of a curve given its curvature 
and torsion) is realized in terms of the solutions of the AKNS system (equivalent to the 
Darboux equations for the natural frame of the curve). Remarkably, the reconstruction of 
the evolving filament requires taking no antiderivatives: given a fundamental matrix solution 
t;X) of the AKNS system, such that $(0,0; A) is a fixed element of SU(2), then the 
skew-hermitian matrix 

1 d$ 

7 (x,t) = $- 1 — (6) 

dA A=0 

satisfies the VFE (JT]), and corresponds to q via the Hasimoto map (26], [29]. (We have identified 
su(2) with M 3 via a fixed isometry, under which the Lie bracket corresponds to -2 times the 
cross product.) Formula (El), known as the Sym-Pohlmeyer reconstruction formula, can also 
be evaluated at a nonzero real eigenvalue A = A . The resulting curve 7 still satisfies ([[]), 
but with a potential that differs from what we started with by the Galilean transformation 
q(x,t) 1 — > e 1 ( ax - a *)g(x — 2at,t), a = — 2A , which preserves solutions of ([3]). Given a closed 
curve of length L, the potential q obtained by the Hasimoto map §Z§ is not necessarily L- 
periodic, but will be related to an L-periodic potential by a Galilean transformation. The 
closed curve may then be recovered using ([HD, but evaluated at A = A . 
The Floquet spectrum of a finite-gap solution. The spectrum associated with an 
L-periodic NLS potential q(x) is defined in terms of the Floquet discriminant 

A(q; A) = Trace($(a; + L, t; A)$(x, t; A)" 1 ), 

the trace of the transfer matrix across one period L, where $ is a fundamental matrix solution 
of the AKNS system. The Floquet spectrum is the set of complex A values for which the 
eigenfunctions of the AKNS system are bounded on the spatial domain: 

a{q) = {A G C I A{q; A) G E, -2 < A < 2} . 

It can be shown that A is a conserved functional of the NLS time evolution, and in fact a 
generating function of the constants of motion. In particular, the spectrum of a given NLS 
potential is invariant under the evolution. 
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Figure 1. The Floquet spectrum of a finite-gap NLS potential. 



In Figure [U we show a typical spectrum for a finite-gap potential, which has a finite number 
of complex bands of continuous spectrum. Those points at which A = ±2 are divided into 
simple points and multiple points, according to their multiplicity as zeros of A 2 — 4. (Finite- 
gap potentials are characterized by having a finite number of simple points.) Multiple points 
are among the critical points for A; for example, double points satisfy dA/dA = 0, but 
d 2 A/dA 2 7^ 0. Real double points are removable [28], i.e., at these A-values the transfer 
matrix has a pair of linearly independent eigenvectors. 

Critical points where A does not achieve its maximum or minimum value are embedded 
within bands of continuous spectrum. It can be easily shown that dA/dA = at transverse 
intersections of bands of continuous spectrum. Those critical points associated with spines 
intersecting the real axis play an important role in the closure condition of the reconstructed 
curve, as well as in our analysis. 

The quasi-momentum differential and the closure condition. At a generic complex 
A, the transfer matrix has a pair of distinct eigenvalues p + (A), p~(A) known as Floquet 
multipliers. The relation between Floquet multipliers and Floquet discriminant is given by 

Thus, p ± (A) are branches of a holomorphic function p which is well defined on a two-sheeted 
Riemann surface S, whose projection ir : £ — > C is branched at the simple points. For 
P G S, we will let 

p(P)=e iLn ^ p \ 

where the function Q\(P) is defined up to adding an integer multiple of 2n/L. Its differential 

dfii = — dlogp, 
lL 

known as the quasimomentum differential, is a well-defined meromorphic differential on E. 
Because dOi changes by minus one from sheet to sheet, each pair of its zeros projects to a 
single A-value, and we will regard its zeros as being located in the complex plane. 

For a given finite-gap NLS potential, the real zeros of dQi (which are the real critical 
points) play a role in the following result by Grinevich and Schmidt |15j : 
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Closure Conditions. A finite-gap VFE solution 7 obtained from the generalized Sym- 
Pohlmeyer reconstruction formula (J5J) at A = Aq is smoothly closed if the reconstruction 
point Ao is (i) a real double point and (ii) a zero of the quasimomentum differential. 

(See also [6] for a derivation of the closure conditions from the explicit formulas for finite-gap 
VFE solutions.) 

Closure conditions must of course be satisfied if one is interested in establishing connections 
between the knot types of closed finite-gap VFE solutions and the spectra of the associated 
NLS potentials. However, such conditions turn out to be difficult to compute, as they involve 
solving implicit equations written in terms of hyperelliptic integrals. 

The central idea of this paper is to approach the problem of constructing closed finite-gap 
solutions of the VFE by starting with an already closed curve of particularly simple type 
(namely a multiply covered circle) and deforming its associated spectrum in such a way that 
both periodicity of the corresponding NLS potential and closure of the curve are preserved. 
By iterating similar isoperiodic deformation steps, we construct a neighborhood of the initial 
curve that consists of closed finite-gap solution of increasing complexity. (Indeed, at each 
step of the deformation, a real double point is "opened up" into a new pair of branch points, 
thus increasing the genus of the Riemann surface by one.) 

A beautiful consequence of this multi-step construction is a labelling scheme that matches 
the deformation data with the knot type of the resulting filament. Our main result, the 
Cabling Theorem (Theorem 15. ip is that every step in this process is, from the topological 
point of view, a cabling operation on the previous curve, with the cabling type encoded in 
which real double points are selected to be deformed into a new pair of branch points. A 
simplified statement of this result is: 

Cabling Theorem. Given relatively prime integers n, m 1; . . . , such that \m,k\ > n > 1, 
let gk = gcd(n, m 1; . . . , m^) for 1 < k < K. Then there exist a sequence of deformations 
(i.e., one-parameter families) of finite- gap potentials q( k >(x; e) of fixed period L, and complex 
numbers A^\e), which are both analytic in e, such that 

(1) qW is of genus k when e ^ 0, q( k \x; 0) = g( fe-1 )(x, e^-i) for some ek-i > when k > 1, 
and q^'(x; 0) is constant; 

(2) the filament (x; e) that is constructed from q^ (x; e) using the Sym-Pohlmeyer formula 
evaluated at (e) is closed of length mr / andis, for e sufficiently small, a (gk-i/gk, mk/gk)- 
cable about 7^(1; Cfe-i)- 

The full statement of Theorem 15.11 explains how these deformations arise from isoperiodic 
deformations of the spectrum, and how the data [n, mi, . . . , uik] determine the selection of 
double points to be opened up. From an argument in the proof of Theorem 15.11 we deduce 
as a bonus that the knot types of the finite-gap filaments so constructed are constant under 
the VFE evolution (see Corollary 15.21) . ultimately confirming and justifying the use of the 
Floquet spectrum as a tool for classifying the knot types of closed curves in an appropriate 
neighborhood of multiply covered circles. (In fact, such curves are approximated by finite-gap 
solutions that are dense in the space of periodic solutions of the VFE [13J.) 

The proof of the main result brings together a variety of tools from the periodic theory 
of integrable systems and the perturbation theory of ordinary differential equations. We 
mention them below while giving a brief description of the paper contents and organization. 
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Section 2 contains the motivation of this work and the framework of our approach. After 
introducing Grinevich and Schmidt's isoperiodic deformation system, we define a special 
type of closure-preserving isoperiodic deformation that reverses "pinching" of the two ends 
of a spine of spectrum into a real double point (homotopic deformations), and show that the 
solution to such a deformation is analytic in the deformation parameter. We then proceed 
to deform off the spectrum of a modulationally unstable plane wave solution (corresponding 
to a multiply-covered circle solution of the VFE), and compute the spatial frequencies of the 
resulting finite-gap solutions. Examples of the curves resulting from successive homotopic 
deformations, and a labelling scheme for their knot types are presented in Section 3. Section 
4 makes use of the completeness of the family of squared eigenf unctions for the AKNS system 
to characterize, to first order, the perturbations of the potential q associated with homotopic 
deformations. The Cabling Theorem and its proof are the contents of Section 5: the proof is 
a combination of explicit perturbative computations involving the squared eigenf unctions, a 
topological argument based on White's formula for self-linking, and a careful analysis of the 
argument of the first order correction to the initial potential, which determines the cabling 
phenomenon and the cable type. The two appendices (Sections 6 and 7) contain a statement 
of the completeness theorem for squared eigenf unctions and related useful results, and a 
proof of the analytic dependence of the potential q on the deformation parameter. 



Inspecting the formula (1361) for a finite-gap NLS solution q(x, t) shows that q is periodic 
in x if and only if the components of the frequency vector V e M. 9 and a real scalar E are 
rationally related. These data are determined by a choice of g + 1 pairs of conjugate branch 
points in the complex plane. Furthermore, because E changes additively when the branch 
points are translated in the real direction, construction of a periodic solution depends on 
being able select the components Vj of the frequency vector. We now describe a scheme for 
deforming the spectrum of a multiply-covered circle which produces arbitrary rational values 
for these components. 

Grinevich and Schmidt [Q] developed a method for deforming the branch points in a way 
that preserves the components of the frequency vector. Such isoperiodic deformations of 
the spectrum were first introduced by Krichever [21] in connection with topological quantum 
field theory, and are naturally related to the Whitham equations in the work of Flaschka, 
Forest and McLaughlin [10]. Although the zeros a>i, . . . , a g+ i of the quasimomentum differ- 
ential dfti are dependent on the branch points Ai, . . . , A2 9 +2, when we incorporate these 
as dependent variables the isoperiodic deformation becomes a system of ordinary differential 
equations with rational right-hand sides: 



Here, ci, . . . , c g+ i are controls, i.e., arbitrary functions of the real deformation parameter £. 
In the case of finite-gap NLS solutions, the Xj and are roots of a real polynomial, and it is 
easily checked that complex conjugacy relationships among these roots (e.g., A 2 j +2 = A2J+1, 
a 2 = «i) «3 = S3) are preserved by these deformations, provided that the controls have 
the same conjugacy relationships as the a^. 



2. Isoperiodic Deformations 





(7) 
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The change in the value of the quasimomentum f2 x at one of the a k under this deformation 
is given by 



Thus, we may preserve the value of Qi(a k ) simply by setting c k to zero. In particular, if 
a k is real and the value of Qi(cek) is such that the Sym-Pohlmeyer reconstruction formula 
(JHJ) yields a closed curve at A = X k , then by choosing the control c k = the isoperiodic 
deformation will produce a closed curve for every £. We will refer to this specialization of 
isoperiodic deformations as a homotopic deformation, since it generates a homotopy through 
the family of smooth maps of the circle into M 3 . 

Instances of homotopic deformation have been observed before. David Singer and the 
second author [17] constructed one-parameter families of closed elastic rod centerlines (which 
correspond to genus one finite-gap NLS solutions under the Hasimoto map) in the form of 
torus knots, terminating in multiply-covered circles at either end of the deformation. We 
can try to generate this deformation using system (JTJ with g = 1. Assume that, say, a 2 
is the real critical point that yields a closed curve. (Necessarily, the other critical point a\ 
must be real.) Then choosing c 2 = and c% = 1 will reproduce part of this homotopic 
deformation. The 'circular' end of the deformation occurs when a\ and one of the complex 
conjugate pairs of branch points (say, Ai and Ai) limit to the same real value as £ decreases 
to some finite time T\. Note that this is a singularity of the isoperiodic deformation system, 
as the right-hand side of (jTJ) blows up as Ai — > a\. In the next subsection, we will examine 
this type of singularity for (JTJ in more detail. 

The genus one solution of (jTJ) discussed above also reaches a singularity in finite forward 
time, when the two a's collide. In our previous paper [B] we showed that this other kind 
of singularity is associated with the elastic rod centerline becoming an Euler elastic curve. 
Presumably, the solution may be continued smoothly through the singularity, although we 
have not investigated this question. 

2.1. Pinch- Type Singularities. We will say that a solution of (JTj) (in arbitrary genus) 
has a pinch-type singularity when exactly two complex conjugate branch points and exactly 
one real critical point a k approach the same real value. (The reason for the name is that 
bringing two branch points together collapses a homotopy cycle on the associated Riemann 
surface, hence pinching one handle on a g-handled torus.) We will limit our attention to the 
case where only the control associated to a k is nonzero. 

Because the system (J7J is invariant under permuting the indices on the branch points, and 
invariant under simultaneously permuting the indices on the critical points and the controls, 
we can without loss of generality let Ai and A 2 = Ai be the colliding branch points, and a.\ 
the critical point, with C\ — 1 being the only nonzero control. Then the system takes the 
form 




dA, 



1 



1 < j < 2g + 2, 



d£ 





1 



2 < k < g+ 1. 



d£ 



at - a k 
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In particular, because A 2 = Ai, 

d(Ai-ai) 1 \ ly^ 2 1 y 1 

d£ 2VAi-ai Y 1 -a 1 ) 2 j-^ Xj - a x ^a k -ai 

showing that, if Ai — a.\ approaches zero as £ \ 0, and the other differences Xj — a\ and 
a k — a,\ have nonzero limits, then lm(Ai — a\) should approach zero like \/2£. 

Proposition 2.1. When the change of variable £ = \t 2 is made in (JSj), the resulting system 
has a solution which is analytic at t = and satisfies 

a k = a° k + 0{t 2 ), l<k<g + l 
Ai = al + it + Oif), 
A2 — Ai, 

Xj = \v + 0(t 2 ), 3<j<2g + 2, 

where a\ is real, and for k > 1 and j > 2, the values a° k , A° are distinct from a\. 
Proof. Let Ai — o>\ = x + iy and z = Re(Ai). Then, because cii is real, 

dz _ ( 1 A. —x 
— = — Re 



d£ \Ai — «i/ x 2 + y 2 

Because the remaining a's are real, and the Xj are in complex conjugate pairs (with A2 = Ai), 
dy 



1 V 1 1 V 1 

d£ \ Ai — a% f--f afc — cti 2 Aj — ai 



-Im 



fc>i * j 

— Im 



2 \Xi — ai Ai — aiJ x 2 + y 



and 



dx 



d£ m [ Ai - ai f^ak-an 2^Xj-a x ) 2 Xj - an ~ a* - <*i " 



Note that the quantity iV is, by assumption, a nonzero analytic function of its arguments at 
their initial values Aj = A°, ot\ = a°, and a k = a k , for j > 2 and > 1. 

We now change to y as independent variable, and introduce the new dependent variables 

x = x/y, z = (z-al)/y, Xj = (Xj - Aj)/y, 5 fc = (a* - a k )/y. 



Of course, if the conclusions of the proposition hold, then these new variables will be analytic 
functions of y that vanish to at least order one when y = 0. In terms of these new dependent 
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and independent variables, the system becomes 

y— = —x + y(l + x 2 )N 
dy 

dz „ „ 
y— = -x-z, 
dy 

dXj ~ y(l + x 2 ) 



xO-aO + y^-z + Xj)' 



da k „ y(l + x 



2\ 



y-r~ = ~ a k 



dy a° k - a\ + y(x - z + a k ) 



This system has a Briot-Bouquet singularity at the origin. Linearizing the right-hand sides 
at x — ~z — Xj — Sfe = gives a lower-triangular coefficient matrix all of whose eigenvalues 
are equal to —1. Thus, there is a unique solution all of whose components vanish when y = 
and which is analytic in y near y = (see, for example, p. 21 in [T8] or Theorem 59 in p3j3) • 
Consequently, x, z, Xj and a k are analytic functions of y which satisfy 

x = 0(y 2 ), z = a° 1 + 0(y 2 ), X, = A° + 0(y 2 ), a k = a° k + 0(y 2 ). (9) 

Because dC,/dy = y(l + x 2 ), it follows that £ is an analytic function of y 2 , and we can 
arrange that £ = when y = 0. Then £ = \y 2 + 0{y A ). Thus, we can define an analytic 
function t of y such that £ = \t 2 and t = y + 0(y 2 ). Clearly, this function is invertible for \y\ 
sufficiently small, so we may replace y by t in (Q, and assert that y is an analytic function 
of t satisfying y = t + 0(t 2 ). Rewriting these equations in terms of the original variables Ai 
and «i completes the proof. □ 

Note that the variable t in Proposition 12.11 will be replaced by the more commonly-used 
deformation parameter e in later sections of this paper. 

2.2. Deforming the Multiply-Covered Circle. Our main idea is to use the isoperiodic 
deformations provided by Proposition 12.11 to create new, higher-genus NLS solutions while 
maintaining the closure of the filament (in effect, reversing the collapse of branch points 
described above in the genus one case). By continuity, the initial value a? used must be 
a real point of the discrete spectrum (hence, a real double point) of the lower-genus NLS 
solution. So, when we carry out several deformations in succession, we must keep track of 
the continuously changing positions of those double points that we want to split later. We 
can do this by incorporating them as extra variables in ([H]), consisting of one critical point 
and a pair of branch points that remain equal as the deformation progresses. (Note that, for 
example, the equality A3 = A4 = «2 is preserved under the deformation ([8]).) 

We will begin the multi-step deformation process by selecting g double points from the 
spectrum of a multiply-covered circle to be the initial values for ax, ... , a g , and begin the 
deformation by setting A2/C = X2k-i = &k for 1 < k < g. (There will be two addition 

1 In the sources cited, the results are given for Briot-Bouquet systems zdyi/dz = Fi(z, y\, . . . , y n ) where 
the linearization has k positive eigenvalues, the resulting solution is a double power series in z and In z that 
converges for \z\ sufficiently small in a sector about the origin in the complex plane, and the solution depends 
on k constants. It is implicit in the statements of these results that, when k = 0, the solution is unique and 
analytic in z. In that case, the existence/uniqueness result can also be proved in a more elementary way 
using the method of majorants. 
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branch points with initial values \2g+i — 1 and \2g+2 — and an additional a g+ i initially 
equal to zero.) The first deformation, which produces a genus one (elastic rod) solution, 
will be continued up to a small value of the deformation parameter. Then, «2 replaces a\ 
as the initial value for the next double point to be split, and so on until all g double points 
have been split. Note that we will assume that the amplitude of each deformation step 
(i.e., the maximum value of £ attained) is sufficiently small so that singularities are avoided. 
(Consequently, all critical points remain real.) As we will see in later sections, it is also 
necessary to assume that the amplitude of each step is small relative to the last one, in order 
to be able to predict the topological type of the resulting filament. 

We will be able to calculate the frequencies for the genus g NLS solution by using the 
fact that the frequencies are unchanged by isoperiodic deformations, and are determined 
by the integrals of certain holomorphic differentials around cycles on the Riemann surface. 
These integrals can, in turn, be calculated by following the sequence of deformations back to 
the multiply-covered circle, and doing residue calculations using the initial positions of the 
double points. 

We begin this calculation with the spectrum of the multiply-covered circle. Using the 
Hasimoto map (T5]), the potential q(x) = 1 corresponds to a circle with curvature k = 2, 
radius 1/2, and length tt. A fundamental matrix solution for the spatial part of the NLS 
linear system 



d</> 

da: 



-iA 
iq 



iA 







(10) 



with q — 1 is given by 



$(x;A) 



COS( LUX 



— sinful) 



- sin (lux) 



- sin [lux} 
cos(c<jx) + — sin(a;x 



to 



The periodic spectrum for the singly- covered circle, computed relative to period L = tt, 
consists of two simple points A = ±i, and infinitely many real double points given by A = 
±vV-l, m = 1, 2, 3, . . .. (In fact, the origin is a point of order four.) For the n-times- 
covered circle we compute the periodic spectrum relative to period L = mr, and this consists 
of the same simple points, but with double points given by 



A = ±a/ {m/n)' 



m 



1,2,3, 



Note that the origin is the only point of the periodic spectrum that can be used as the value 
A producing a closed curve via the Sym-Pohlmeyer reconstruction formula Thus, zero 
will be one of the initial values for the critical points in the first step of the deformation 
process, and the corresponding control will always be set to zero in order to maintain closure; 
we will use Ao to denote this 'reconstruction point', which will move along the real axis during 
the deformation steps. 

Let (3\ < fa < ■ ■ ■ < P g denote the choice of double points of the spectrum of the n-times 
covered circle to be opened up. (For the following calculation, it is not relevant in what order 
they are opened up.) We will now determine the relationship between these double points 
and the periods of the genus g solution obtained once all g points have been opened up. 
According to the construction for finite-gap solutions set forth in [2] and [5], the frequency 
vector [Vi, . . . , V g ] is Airi times the first column of the inverse of the matrix A defined by 



Ai 



a k 
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where ( is related to A by the defining equation of the hyperelliptic curve, 

2g+2 

c 2 = n( A -^)> 

3=1 

and each cycle circles around (in a clockwise fashion, on the upper sheet!) the pair of 
complex conjugate branch points into which /3 k has been split (see Figure ED- 

First, consider the case when all /?'s are negative. Then, as we run the deformation steps 
backward, the cycles Ol, . . . , CLg become loops around (3i,...,(3 g , and the only remaining 
branch points are ±i. The denominator £ in the integrand limits to — \/A 2 + 1 nf=i(^ — Pj) 
along the a-cycles, taking the square root as positive along the real axis. A residue calculation 
then gives 



lim/r^. 



27Ti/?f J 



<M+iU.(Pk-Pi) 

(Note that the minus sign in ( is offset by the clockwise orientation of the a-cycles.) We 
may factor this matrix as 



N = 2m 



'PI 1 

PI 2 



PI 1 ' 

pr 



l/D 1 








l/D 2 









Di 



pt + ^H^-fo 



(12) 



Its inverse is 



1 
2^i 





. 


. 







D 2 . 


. 


1 













. 


• D 9 



lj<k;j,k^l(Pj Pk) 

(-1) Il j<k;j , k MPi ~ A) 

(-i) 2 n i<W3 (^--^) 



(Only the first column is needed for the matrix on the right.) Thus, 

2, 



Pj + 1, 



l<3<9- 



(13) 



Next, suppose that /?i < . . . < (3k < < (3k+i < ■ ■ ■ < P g - Because the branch cut 
between i and — i, the roles of upper sheet and lower sheet are switched an extra time, and 
the residue calculation gives 

N k} 1 < k < K, 

lim A k 



K + 1 < k < g, 



where N k is the fcth column of the matrix defined by (TT2"]) and A k is the fcth column of matrix 
A. Solving for lim A -1 gives the following frequency formulas: 



Yi 




1<J<K, 
K+l<j<g 



(14) 



2 Along the real A axis, we initially label as the upper sheet that containing the point oo+, where A 9+1 /£ 
tends to +1 as A — > +oo, but as one proceeds from right to left along the real axis in the A-plane, the roles 
of upper sheet and lower sheet are exchanged along branch cuts that run between each branch point and its 
conjugate, parallel to the imaginary axis. 
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3. Examples of cabling operations 

Comparing the frequency formulas (fTBl and (1141) with formula (lllj) for the spectrum of the 
multiply-covered circle shows that when the double points (3j are chosen as members of that 
spectrum, the frequencies attained by the deformations have rational values. Accordingly, 
we will introduce a notation for these deformations that allows these values to be read off 
easily 

We will use the notation 

[n;mi, . . . ,m g ], \rrij\ > n 

to indicate the result of g successive homotopic deformations of the n-times covered circle, 
opening up the real double points whose starting position is (3j = — sign(mj) (rnj / n) 2 — 1 
in the order in which the rrij appear in the square bracket sE (The minus sign is incorporated 
here to compensate for the sign in f)14p .) As mentioned before, we will assume that the 
amplitude of each deformation step is small relative to that of the previous step. We will 
also assume that the numbers in square brackets rrii are relatively prime, so that the double 
points selected do not all belong to the spectrum of a multiply-covered circle for a smaller 
value of n. In this notation, it is easy to see that [n; m] gives frequency V = m/n, [n; mi, 
gives frequency vector [mi/n, m2/n], and so on. These frequencies determine the length of 
the corresponding filament as follows: 

The finite-gap NLS solution q(x,t) takes the form 

q(x, t) = Aexp(-iEx + iNt)6(iVx + iWt - D + r) /9{i\x + iWt - D), 

where A,E,N are constants, V is the frequency vector, D and r are constant vectors in 
C 9 , and 9 is a Riemann theta function with periods that are 27ri times an integer vector. 
(For more details, see §0) The phase factor exp(— iEx) may be removed by an appropriate 
gauge transformation of NLS, and then, assuming the components of V are rational, the 
period of the potential q is the least common (integer) multiple of the periods 27r/Vj. The 
Baker- Akheizer functions, quadratic products of which are used to reconstruct the Frenet 
frame of the filament, take the form 

rj,. = d(P) expiixtt^P) + itn 2 (P))^(iVa; + iWt + Di(P))/0(iVa; + iWt - D), z = 1,2, 

where P a point on the Riemann surface lying over the reconstruction point Ao, fii,2 are 
certain Abelian integrals on S, and C\ t 2 and D12 are constants that depend on P. Under the 
homotopic deformation process we have described above, the value Qi(P) is preserved, so 
can be obtained by calculating its limit as we run the deformation process backwards to the 
multiply-covered circle. Again, we first calculate this assuming that all the limiting values 
Pj of the double points are negative. Then, because 

9+1 

dfii = ^— dA, 

the limiting value of this differential, on the upper sheet above the origin, is 

lim dCli 



3 However, it matters in which order the deformations are done, for it is easily checked that if Vj is the 
vector field on C 3s+3 defined by setting Cj — 1 and all other controls zero in the system ([7]), then the vector 
fields Vj and v^ do not commute when j 7^ k. 
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Integrating from the limit — i of the basepoint for (which is by convention the point in 
the lower half-plane not enclosed by the basis a-cycles) to the origin, which is the limit of 
the reconstruction point, gives fli(P) = 1. Next, assuming that the least positive double 
point is Pk+i, 

lim dfi! = (-1)^_SL. 

Integrating this from the basepoint to the origin gives fl\(P) = (— 1) 9 ~ K . Thus, the factor 
exp(2ix^i(P)), which occurs in quadratic products of Baker functions, has period 7r. Since 
this factor is multiplied the theta functions, the period of the Frenet frame is the least 
common multiple of tt and the numbers Zn/Vj for 1 < j < g. Because the reconstruction 
point is chosen so that the integral of the Frenet unit tangent vector is zero, this period is 
also the length of the filament. 

As we will eventually show, the selection of frequencies determines the knot type of the 
resulting filament as an iterated cable knot. This can be viewed as a generalization of the 
work of Keener [TJ] , who showed that if one adds to a circle of length L a small perturbation 
of period (n/m)L, and a closed curve results, then the perturbed circle is a (n, m) torus knot 
(i.e., it covers the original circle n times lengthwise, and wraps around the circle m times). 
The Figures [2] through H] show examples of our iterated cable construction. 

4. Isoperiodic Deformations and Squared Eigenfunctions 

In this section we show how, for an NLS solution go of a given period, the only isoperiodic 
deformations that preserve the reality condition r = — q (see $6] for notation) and open up 
just one real double point Ao, while leaving all other points of the discrete spectrum and the 
critical points unchanged at first order, must be linear combinations with real coefficients 
of (<^ + )i + ( ( P2 + ) 2 and i[(^ + )i — { i p2 + ) 2 }i written in terms of the components of a Bloch 
eigenfunction cp + evaluated at Ao- The proof requires computing the first order variations 
of the discrete spectrum (simple and multiple points) and of the critical points. For the 
purpose of this paper, it is sufficient to discuss the computation of the first order variation 
of a real double point; the cases of simple and critical points can be treated analogously. 

Assume that (go + egi, + e r i) is a periodic perturbation of a NLS potential qo = (go, To) 
of a fixed period L, and assume that the pair (</> , Ao) solves the AKNS eigenvalue problem 
at qo, where Ao is a real double point in the spectrum of the unperturbed potential qo. 
Let 4> = 4> + e<p 1 and A = Ao + eAi be the corresponding variations of eigenfunction and 
eigenvalue. At first order in e, we obtain 

- {-\rl it) < 18 > 

where C the spatial operator of the unperturbed AKNS linear 



dx \-ir i\ 0/ 
system at A . 

Together with the homogeneous linear system Ci(f> = 0, we consider its formal adjoint 

with respect to the L 2 -inner product (u, v) = / u • vds: 

Jo 

d f-iXo -if 



C?il> = 0, with £f - , . _ 

dx V igo 1A0 
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Figure 2. Results of the two-step deformation [4;— 10,7]. The end of the 
first step is the left-hand (2, —5) torus knot shown at top, and the result of 
the second step is the knot shown at bottom, which is a right-hand (2, 7) cable 
on the torus knot. Approximate diagrams of the Floquet spectrum appear at 
right; in these, the reconstruction point is marked by a diamond, while double 
points are not shown. 



(The "formal adjoint" is computed neglecting boundary conditions. It turns out that the 
solvability condition of system (Tl5]) is correctly expressed in terms of the formal adjoint, 
as discussed below.) Taking the inner product of both sides of (|15p with a solution of 
CfiJ) = leads to the solvability condition for system ffTBT) . requiring that its right-hand side 
is orthogonal to any solution of the homogeneous formal adjoint system. 

Since A is assumed to be a real double point, it is removable [28]. This means that the 
space of solutions of C±(j> = is spanned by the two linearly independent Bloch eigenfunctions 



14 



A. CALINI AND T. IVEY 




Figure 3. Results of the two-step deformations [4; —6, —13] (in the middle) 
and [6; —9, 10] (at bottom). In both cases, the result of the first step of the 
deformation gives the left-hand trefoil shown at top. The [4; —6, —13] is a left- 
hand (2, —13) cable on the trefoil, while the [6; —9, 10] is a right-hand (3, 5) 
cable on the trefoil. 
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Figure 4. Result^of the three-step deformation [8; -12, 10, -17]. The first 
step gives the left-hand trefoil shown at top, the second step (shown in the 
middle) gives a right-hand (2, 5) cable on the trefoil, and the third step (shown 
at bottom) gives a left-hand (2, —17) cable on the previous cable. 
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evaluated at \q. It is easy to show that a basis for the solution space of the adjoint system 
Cfi/) = is given by 



Jcp- 



¥>2_ 





Jif 



where J is the symplectic matrix I ^ ^ 1. Substituting O = a(p + + bip , a general solution 

of Ci4> = 0, in system ( Fl5l) . we compute the solvability condition to be the following system 
in the unknowns a and b: 



2iAi J ipf(p% dx + i(^qi, J A \ a+ iXi J {<pftp 2 + <ft<Pi ) da: + i ^qi, J 



< PI<P1 



a + 


2i\ 1 [ 




Jo 



-^2^2 , 



6 = 0, 



6 = 0, 



where the various components of the Bloch eigenfunctions are evaluated at the double points 
A . 

For a removable double point A^, the following identities hold 

-L r L 



WW I a. 



dx = 0, 



WW| Ad 



dx = 0, 



as shown, for example, in [15J. Thus, the condition for existence of a nontrivial solution 
(a, b) of the above system (amounting to the vanishing of the determinant of the associated 
matrix) provides the selection of the perturbation of the affected double point. A simple 
calculation gives 





( ift^Pl \ 






( (^) 2 ' 


) 


J( ql ' J 


( (^) 2 \, 




Jo (<Pl<P2 +<pt<Pl)\ 


, dx 



Ai 



where qi denotes the vector (gi,ri) T . 

We now recall that certain vectors of quadratic products of components of the Bloch 
eigenfunctions (the squared eigenfunctions) form a basis for Lp er ([0, L], C). (See §6] and 
references therein.) Moreover, we can write a generic periodic perturbation q! as a linear 
combination of the elements of the basis of squared eigenfunctions: 



° s 

simple points \ s 



+ 



E 



s critical non-per. pts A c 



+ 



E 



remov. double pts. 



<Pl<Pl 

(<pp 2 



2 i + e 



(v?r) 2 N 

1 1 -tea )\ 



(We have assumed for simplicity that the potential qo possesses no non-removable double 
points and no periodic points of multiplicity higher than two.) 
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Because of the biorthogonality property of the squared eigenfunctions (see §|6]), we compute 




with normalization coefficient (see Lemma 16.31) 




JVA(Ao)A"(A )|WV 



Ac 



2i 



I 2 
I Ao 



where VF[<£ + , cp~~] denotes the Wronskian of the two Bloch eigenfunctions. Because the 
normalization coefficient is non zero, it follows that the real double point Ao will split at first 
order and none of the remaining double point will if and only if the perturbation qi contains 
only terms of the form 



K) 2 



+ e r " 



and none of the terms e \ 



contains terms of the form 



K) 2 

-(v4) 2 



Ao 



-(<P~2 



Ao 



or 



respectively. Thus, because we are 



to? 
-to) 2 

Analogous results can be deduced by computing expressions for the first order variation 
of simple and critical points: they will move at first order if and only if the potential qi 

( P\ ( Px_ 

assuming that these points do not move to first order, then these terms do not occur in qi. 

Since Ao is assumed to be real, the Bloch eigenfunctions possess the additional symmetry 
<p~ = Jip + , thus a perturbation q x that splits only A , leaving the rest of the discrete 
spectrum and critical points invariant at first order must be of the form: 

"• 2 + e -(<^ 2 



Ao 



Finally, requiring that the focussing reality constraints r\ - 
condition 

(4 -0(v£) a - (4 - e )K) 2 = o, 

and, given the linear independence of ((fit) 2 and (<^2~) 2 over ^> one obtains 
some complex constant c. We summarize these results in the following 



-qi be satisfied leads to the 



c for 



Proposition 4.1. Suppose the periodic potential go undergoes a smooth perturbation q = 
% + e Qi an d its spectrum deforms isoperiodically in such a way a unique real double point X 
splits at first order, while all other points the discrete spectrum as well as the critical points 
remain unchanged up to first order. If, in addition, the deformation preserves the reality 
of the potential, then the perturbation must have the form qi = Re (c)[(^) 2 + (y?2~) 2 ]| A( + 
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i"m(c)[K) 2 -(^ 



\2l 



2 



in the following pair of simple points 



, with c G C. Moreover, the affected double point splits at first order 

Ao 



il=J , ±i£ | cf wtrf#H +0(£!) . 

/„(lrfl 2 -|rfl 2 )dz 

5. The Cabling Theorem 

In this section we will assume that g is an NLS potential of period nir, obtained by a 
sequence of g homotopic deformations from that of the n-times covered circle, notated as 
[n; mi,..., m g ), and that Ao is a double point of the jri7r-periodic spectrum of go whose orig- 
inal position (in the spectrum of the circle, before the deformations) was —a/ (m/ (jn)) 2 — 1 
times the sign of m, where m is an integer whose magnitude is greater that j and which is 
relatively prime to j. We will show that, if we perturb go as m Proposition 14.11 then the 
perturbed curve 7 = 7 + e7i + 0(e 2 ) is an (j, m)-cable on 7 for sufficiently small e. 

We also need to assume that g is sufficiently close to the plane wave potential; we will be 
more specific about this assumption at the end of this section. 

5.1. Perturbed Potentials and Perturbed Curves via the Sym Formula. We will 

begin by calculating the perturbation of the curve. We assume that 

q(x,t;e) = q + eq 1 + 0(e 2 ) 

is a one-parameter family of NLS solutions, and 

$(x,t;A,e) = $ + e$i + O(e 2 ), 

is a fundamental matrix solution for the AKNS system at (g, A). To be specific, we will 
assume that $(0, 0; A, e) is the identity matrix. Moreover, when A is real, we can assume that 
$ takes the form [<fi, J<p\. We will also assume that gi is expressed as a linear combination 
of squared eigenfunctions for the AKNS system at (g , A ), where A is a double point for g . 
(This assumption was justified, for isoperiodic deformations that open up a double point, in 
the previous section.) 

For use in the Sym-Pohlmeyer formula ([H]), we will need to calculate $i at A = A , 
the reconstruction point. (Although this reconstruction point will vary under isoperiodic 
deformation, Proposition 12.11 ensures that it is fixed up to 0(e 2 ).) For the sake of brevity, 
we will use the abbreviations 

<Pi = <Pi(.Qo, Ao), = v4(<7o, A ), 

where ip + is one of the Bloch eigenfunctions, as normalized in (1521) . Then the results of §H 
imply that 

g 1 = c(^ 1 ) 2 + c(^ 2 ) 2 (16) 
for some complex constant c. We will also let 

-<*■*■•*>=(£ -£)• < 17 > 

where we will leave A arbitrary (but real) for the moment, but later set A = A . 

Setting Bi = $x($ ) 1 an d taking the e term on each side of the AKNS system gives 
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(see j5] for more details). Expanding the inverse, we get 



d.r 



D 



1 ( 410102 + 410102 



2 2 

qx<j>x - 4i0 2 



4i01 - 4i0| 



-qn 



h ~ 510102 



(18) 



where D = det $o = |0i| 2 + 102 1 2 , which is independent of x and t. 

To compute Bi, we need to take antiderivatives of the entries of the matrix appearing 
on the right-hand side of ([TBI . (Because the matrix takes value in su(2), it suffices to find 
antiderivatives for entries in the first column.) To do this we will use special properties of 
solutions of the squared eigenfunction system fT28|) . ff29l) Recall that if <p,4> are two vector 
solutions of the AKNS system at the same A-value, then the construction 



0® 4> 



|(0102 + 0102) 
0101 
-0202 



(19) 



gives a solution of the squared eigenfunction system (|28|) . (|29|) (see §6.11 in the Appendix). 

Furthermore, if (/, g, h) and (/, g, h) are solutions of this system at (go, A) and (go, A) respec- 
tively, then 



— {gh + gh- 2ff) = 2i(A - \){gh - gh). 
ax 



(This can be easily verified from the system of ODE's 
example, the proof of Theorem 2.1 in [21].) 

To get the top left entry in (fl8|) . let 'g — qx, h — — q 1 , A = A , and g 



(20) 

satisfied by (f,g,h); see, for 

0102, h ^ 



2Vx 



(the latter arising from the construction </> (g) (J(f>), giving / 

1 



/ 



410102 + <?10102 dx 



410102 - <7l010 2 



2i(A - A 

where, using the formula ( TlBT) for qx, we have 



KI0i 12 
(I0il 2 - 



b 2 \ 2 )). Then 



To get the bottom left entry in fflBl (up to a factor of minus one), we keep f, g, h the same, 



and change to g = <p\ and h = — 0| (arising from the construction cf> ® <f>, giving / 

L2 -= ±2 



Then 



4i0 2 - 4i0i dx 



1 



2i(A - A ) 

So, up to an additive constant, we obtain 



4l01 + 4102 - 20102/) • 



Bx 



1 



2(A-Ao)D 



-410102 + 410102 - (|0l| 2 -J_02| 2 )/ 
-4l0? - 4102 + 20102/ 



_ 4101 +4l02 + 201 2 / 
410102 - 410102 + (|0l| 2 - I02I 2 )/, 



Substituting $ = §o(I + eBi + 0(e 2 )) into the Sym-Pohlmeyer formula gives 



7 = 7o + e7i + 0(e 2 ), 7o = $, 



dA ' 



7i 



dfii 



+ ho,Bx\. 
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(Recall that we are identifying matrices in su(2) with vectors in M 3 .) But B\ = f $ 1 Qi&o cLc, 
and Qi is independent of A, so 



dA 



= y [(Si^To] da: = [£?i,7o] " y dx, 



where integration by parts is used in the second line. Therefore, ~y 1 = J[T, B\\ dx, where 
r _ r i/"i OA _ 1 f-i(|0i| 2 - |0 2 | 2 ) 2i0,0 2 A 

We now set A = A , the reconstruction point that generates the closed curve 7 from potential 
q . Thus, because J T dx = 0, the additive constant in B x does not change the value of j 1 . 

Note that the coefficient of / in Bi is an exact multiple of the matrix T. So, 



2D 2 (A - A ) 



<h\ 2 ~ l^ll 2 2 0102 ^ (-q\<t>\<t>2 + 910102 9101 + 9l02 

20102 |0l| 2 - |0 2 |7 ' V -9l0 2 -gi02 910102-^10102 



1 ( -910102 - 910102 -910J + 9102 



D(A - A ) V 9i0? - 9i02 -9i0i02 - 9i0i<^2/ 
Amazingly, the quantities we have to integrate are the same as before; in fact, 

-1 dSi 



[T, Bi] = 



(A - A ) da; ' 



Therefore, 

— AO 



.Zl A, /0102 0M ^77^0102 02 V f^l02| 2 -|0l| 2 ) 2<M 



2 



2(A - \ fD \ qi V-0 2 -0 X 2 J + Ql V -0? 2 J 2 J + ; V 20!0 2 |0i| 2 - |0~' 2 
When we expand this in terms of the A -natural frame of 7 , comprising T with 

then we obtain the first order perturbation term in the expression for the curve 

71 = 2(A - 1 Ao) 2 (~^ T + Re(gi) Ul + ' m(gi) f/2 ) ' (21) 

5.2. Natural Frames Unlinked. Now that we have a nice expression for the components of 
the perturbed curve in terms of a natural frame, we will show that under some circumstances 
the perturbed curve 7 + ey 1 forms a cable around the unperturbed curve 7 . To determine 
the type of the cable correctly, we need to know if the natural frame itself winds around 7 . 
We will show that, for sufficiently small e, the curve 7 + eU\ is unlinked with 7 — provided 
that 7 has self-linking number zero and has no inflection points. Note that, because the 
self-linking number is a discrete invariant — either an integer or a half-integer — it does not 
change under deformation, and is multiplied by the integer n when we pass from a curve 
to its n-fold cover. Therefore, any curve that is obtained from the circle by a succession of 
multiple-coverings and deformations will also have zero self-linking number. 
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According to White's formula [31], the linking number of 7 and 7 + eU\ 
Lk( 7o , 7o + el/i) = Wr( 7o ) + ^ J {T x ^)d^. 
Because dt/i/dx = — fciT + A C/ 2 (where a; is arclength) and U 2 = T x f/i, then 

Lk(7o,7o + ^i)=Wr(7 ) + ^, 

where L is the length of 7 . Meanwhile, the writhe Wr(7 ) is related to the self-linking 
number by Pohl's formula [25J 

SL( 7o ) = Wr( 7o ) + ^ J rdx. 

Taking SL(7 ) to be zero, we get 

Lk( 7o , 7 + eUi) = J (A - r)dx. 

We can relate the last integrand to the unperturbed potential go = §(&i + 1^2) written 
in terms of the natural curvatures, by making use of the Frenet equations and the natural 
frame equations. Differentiating each side of the expression 

K N = jfejt/i + k 2 U 2 

for the Frenet normal in terms of the natural frame vectors, and cancelling out tangent terms 
on both sides gives 

k x N + ktB = (h) x Ui + {k 2 ) x U 2 + hA U 2 - k 2 A Ui. 
Next, dotting both sides of the with kB = T x k,N = k\U 2 — k 2 U\ gives 

K 2 T = k 1 (k 2 ) x - k 2 (k x ) x + Kqk 2 . 

So, 

« fki(k 2 ) x - k 2 (k 1 ) x \ 

Ao-r = - ( — 2 I = -lm((g )*/gb). 

Now, the value A is chosen so that q is periodic along the curve. So, provided that q is 
nonvanishing along the curve (and, this is true if 7 is sufficiently close to a multiply-covered 
circle) then the integral is zero, and the natural frame is unlinked. 

Thus, the knotting of the perturbed curve about 7 is completely determined by the 
behavior of q±. 

5.3. Monotonicity of Argument. In Section 4 we showed that if q\ only contains terms 
associated to the double point being opened up at first order, then it must be of the form 
( I16p . We will establish that the argument of this function is monotone in x by calculating it 
for the multiply-covered circle, and invoking continuous dependence on the potential. 

As an exercise, one can calculate the following Baker eigenfunction for the plane wave 
solution q = exp(2it) for real values of A: 

exp(it) 



, . , exp(ikx + 2ik\t) 
*l>(x,t;\) = 



(k + A) exp(-it) 



k + X + 1 

where k = VA 2 + 1. This expression extends to the spectral curve fj 2 = A 2 + 1 associated 
to the plane wave solution by replacing k with /a, and coincides with the expression of the 
Baker- Akhiezer eigenfunction normalized as in [6]. 
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We then calculate that 



cipl + cijj 2 = exp(2it)- 



(k + X + 1) 2 

where u := cexp(2iA;x + 4iA;At). (Note that A and k are assumed to be real.) Because 
du/dx = 2iku, then 



log (cip 2 + cip\ J = 2ik 



(1 + (A; + A) 4 )|«| 


2 + (A; + A) 2 (« 2 -?Z 2 ) 




u + (k + A) 2 w| 


2 



dx & V ri r V M + (fc + A) 2 u' 

Then 

lm ^,o g (^ + -^)=2,Re^|±|| = 2*Re 

= aM . (lzi*±m^3 = _ 4 a|cpM±|±^!), ,22) 
|w + (A; + A^-up I it + (A; + A)nt| 2 

where in the last step we use the fact that k 2 = 1 + A 2 to write 1 — (k + A) 2 = — 2A(/c + A). 

The sign of this quantity explains why opening up double points where A > leads to 
left-handed cables, and A < leads to right-handed cables. In the next section, we will 
determine the type of this cable. 

5.4. Cable Type. Given a knot K\ : S 1 — > M 3 whose image lies inside a solid torus T 3 , 
and a knot K 2 : S* 2 — > M 3 , we can form a satellite knot on K3 = t o K\, where r is a 
diffeomorphism from T 3 to a tubular neighborhood of the image of K 2 . In particular, when 
K\ is a (p, g) torus knot, then K 3 is a (p, g) cable on K 2 . 

Because the natural frames of 7 are unlinked, we can use them to define a diffeomorphism 
5 from S 1 x D 2 to a tubular neighborhood of 7 , given by 

5 : (z,r,0) i-> 7 (^) + (rcos6)Ui(z) + (r sin 6) U 2 (z) , 

where 2, r, # are cylindrical coordinates, with r less than one half of the minimum radius 
of curvature of 7 , and we think of S 1 x D 2 as the quotient of the solid cylinder by the 
equivalence relation z ~ z + L. 
By ( jUD , the perturbed curve is 

7(0;) = 7o (a;) - 2(Ao € _ ^ )2 (-i/T + Re( gi ) fA + lm( gi ) C/ 2 ) + 0(e 2 ). 

Consider the modification 

7(2) = 7 (x) - — - 6 (Re(gi) C/i + Imfo) C/ 2 ) , 

which, by continuity, will have the same knot type as 7(x) for e sufficiently small. Because 
qi is composed of quadratic products of eigenf unctions, which have the form (fTUI) . the period 
of qi is the lowest common multiple of L = nn and 7r/f^ 1 (A ). Because the isoperiodic 
deformations preserve the value of ^i(Ao), we can compute it by running the deformations 
backwards to the multiply-covered circle. Specifically, suppose that [n; mi, . . . , m g ] is the 
deformation scheme that produced go- Then, upon running the deformation steps backwards, 
the basepoint for fii limits to — i. So, 



m 



'V (m/0'")) 2 -l \ 

n 1 (A )=/ ^f^ dA 

VA 2 + 1 



Hence, 7r/Oj(Ao) = jL/\m\, and the period of gi is jL. (Recall that j and m are coprime.) 
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The image of 7(2) under the inverse of 8 is the curve 

Z = X} r = e|g 1 (x)|/(2(A - A ) 2 ), 6 = aig(- qi (x)). 

Note that r is never zero along this curve, so that it never crosses the z-axis. Consider the 
map 9 : S 1 —>■ S 1 defined by x 1— > arg(gi(x)), where the domain S 1 is R modulo jL and the 
codomain S 1 is K modulo 2n. For a fixed c, 

gi = ct/; 2 + cipl = c (ipf(q , A )) 2 + c(^(g , A )) 2 

depends continuously on the location of A and the potential q . Thus, the degree of the 
mapping 9 is unchanged under the multi-step isoperiodic deformation process. We can 
calculate its degree by substituting k = —m/(Jn) and A = — sign(m) i/F+1 into formula 
( 122]) . Integrating that formula shows that 

^^-H^uM)- (23) 

So, the period of arg(gi) is 7r/\k\ = jmr/\m\. Therefore, the degree of mapping 9 is m, taking 
the minus sign in (123]) into account. 

Because the curve 8~ l ('y(x)) winds around the z-axis \m\ times (in a counterclockwise 
direction when m > and clockwise when m < 0) we see that it is a (j, m) torus knot in 
S 1 x D 2 . Thus, the curve 7(2) is a (j, m) cable about 7 . 

5.5. Conclusions. We may now state our main result in full: 

Theorem 5.1. Let n, m\, . . .rriK be relatively prime with \rrik\ > n > 0, and let gk = 
gcd(n, mi, . . . , rrik). Let (*) fe be the isoperiodic deformation system with variables oti, Xj for 
< £ < K , 1 < j < 2K + 2, controls = 1 ; q = for £ ^ k, and with change of variable 
i=\t 2 . Then: 

1. For each k between 1 and K there exists pk > such that (*)& has analytic solution 
{af \ Xj) for |e| < pk satisfying, when k = 1, 

a?>(0) = 0, A«(0)=i, A«(0) = -i,^(0) = -^(0) = 1 , 

«f } (0) = A$ +1 (0) = A$ +2 (0) = - signK)V(m,/n) 2 -l, 1 < £ < K. 
and when k > 1 fmft depending on the choice of e k -\ G (0, Pk-ij) 

af >(0) = a^VO, A? } (0) = A^Vi), ^f±i(0) = -^(0) = i 

For eac/i t/iere exist finite- gap potentials q( k >(x; e) which are nil-periodic in x, analytic 
in e, and for which the simple points are A$ fe) (e) ; 1 < j < 2k + 2. Then the filament 7(^(1; e) , 
constructed from using the Sym-Pohlmeyer formula at Ao = ckq (e) ; is closed of length 
nn/gk and is a (gk-i/gk,mk/gk) -cable about 7"(i; €&) . 

The argument given in Section [7] implies that for each k there is a deformation q^ k \x, t; e) 
of NLS solutions whose spectrum matches the deformation of the branch points at any time. 
The above theorem, applied at any fixed time to, implies that the corresponding filament 
7(^(2, t ; £jfc) has the desired cable type. (Note that the plane wave potential q(°> at time t 
only differs from that at time zero by multiplication by a unit modulus constant.) Thus, we 
have the 
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Corollary 5.2. The knot type of (x, t ; e^) is fixed for all time. 

The local nature of the VFE can, in general, cause the knot type to change in time. Indeed, 
it is relatively easy to construct solutions with changing knot type by taking Backlund 
transformations of genus one finite-gap solutions. Thus the Corollary implies that we can 
nonetheless construct a neighborhood of the multiply-covered circle, within the class of finite- 
gap VFE solutions, which consists of filaments whose knot is preserved. 



6. Appendix: Completeness of the Squared Eigenfunctions 



In this Appendix we summarize the main properties of squared eigenfunctions for the NLS 
equation, in particular their connection to the linearized NLS equation, their biorthogonality 
property, and the L 2 -completeness of a suitably chosen periodic subfamily. 

We rewrite the NLS equation as a dynamical system 

(24) 



ir t 



-q xx + 2q r, 
f xx 2r q, 

for the pair (q, r) in the phase space V = Hp er ([0, L), C 2 ) C L 2 ([0, 1], C 2 ) of periodic, square 
integrable, vector-valued functions of x, with square integrable first derivative. The inner 
product is taken to be the one of the ambient space: 

(f , S) = j J\fi{x)9i{x) + h{x)W)\ to, (25) 

where f = (A, /a), g = (fll, 92) are in V . 

The focusing/defocusing reality condition in achieved by restricting to one of the invariant 



subspaces V ± 
inner product 



{(q, r) G V\r = ±q} on which the inner product ff25l) becomes the real 



u(x)v(x) dx 



where f = (u, ±m), g = (v, ±v) are elements of V^. 

For a given NLS potential q = (q, r) G V, we consider the linearization of the NLS system 
(|24l) . obtained by replacing q ^ q + u, r — > r + v and retaining terms up to first order in u 
and v, 



hit + u xx — 2q v — Aqru = 
iv t — v xx + 2r 2 u + Aqrv = 0. 



(26) 



For fixed time t , we regard the pair (u, v) as an element of the ambient space £ 2 ([0, L], C 2 ) 
endowed with Hermitian inner product given by (1251) . 



6.1. Squared Eigenfunctions. If <fi and <fi are solutions of the NLS linear system fl33l) at 
the same value of A, then the triplet 



0® (j> 



>102 + 0102) 

4>x4>x 

-0202 









9 




h 



(27) 



solves the squared eigenfunction systems: 
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7" 
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h 















2(2iAg - 
2(2iAr + r x ) 



— ir — ig 
2ig -2iA 
2ir 2iA 

-(2i\r + r x ) -(2iAg-g a 
-2(2iA 2 + igr) 

2(2iA 2 + igr) 



{21 



(29) 



solves the linearized 




Proposition 6.1. (1) (linearization) The pair (g,h) = (0i0i, 
NLS equations (|26|) . 

(2) (biorthogonality) Suppose that A 7^ are distinct points of the periodic/ antiperiodic spec- 



trum of the spatial linear operator L\ 



d 

icr 3 — + 

dx 



g 
r 



Z-et (0(A), h(X)) and (g(fj,), h(ji)) 



be solutions of the squared eigenfunction systems. Then, 




0. 



(30) 



6.2. A basis of squared eigenfunctions. Given a finite-genus NLS potential q, one can 
construct a periodic subfamily of solutions of the linearized NLS equation which, for fixed 
t, is also a basis for L 2 ([0, L], C 2 ). Its elements are squared eigenfunctions constructed from 
Bloch eigenfunctions evaluated at a countable number of points of the spectrum of q. 

The Bloch eigenfunctions are common solutions of the AKNS system and the shift oper- 
ator, i.e., they satisfy the additional property 



<p ± {x + L)=p ± {\)^ 



where p ± (A) are the Floquet multipliers 



A(A) ± VA 2 (A) - 4 



(31) 



A useful normalization is obtained by selecting (see 

[M 12 (A)Y 



e 4 



e 4 



) 

(P + (A) 



M n (A))Y 2 ] , 
[M 12 (A)Y a + (p-(A) - M n (A))Y 2 ] . 



(32) 



p-(A)-M n (A) 
M 12 (A)(A 2 (A)-4) 

p+(A)-M n (A) 
M 12 (A)(A 2 (A)-4) 

so as to satisfy the useful symmetry 

ip~{x;X) 

In the expressions above, M^-(A) are the entries of the transfer matrix M(x, t; A) = $(x + 
L, t; A)$(x, t; A)^ 1 across a spatial period L, and the vectors Y 4 are the columns of the 
fundamental matrix solution of the AKNS system fl2"B"|) . 

For simplicity, we consider an NLS potential (q, r) the critical points of which are of alge- 
braic and geometric multiplicity two. (Higher order critical points introduce some technical 
difficulty, while critical points of geometric multiplicity one can be dealt with in a simple 
way.) We define the following family of squared eigenfunctions. 



J<p+(x;\). 
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(A) At the real and complex double points {A^}: 

[ } ~{ -(v4) 2 (*; A) J x _ xd ' " V -fo" ) 2 (*; A) j v v , 

i i 

(B) At the finite number iV of critical points {A^} that are not elements of the peri- 
odic/antiperiodic spectrum of q: 



f( c . fc ) 



-(<pt<P2)(x;X) 



(C) At A" (appropriately selected) simple points {A*} of the spectrum (half the number of 
branch points of the associated Riemann surface): 

Then 

Theorem 6.2 (Completeness). The family 

{f ( +J ) if <-,i, } U{f <^» } y {fM} 

is a basis of L 2 ([0, L], C 2 ). 

A proof of this fact can be found in the preprint [7j, following the roadmap developed 
in [9] for the sine-Gordon equation and incorporating ideas from [20]. We also mention the 
following result: 



Lemma 6.3 (Biorthogonal Pairing). The families {f( +J )} and {Jf(~' fc )}; {f( '^}\j\>N an d 
{ Jf(+' fc )} ; are biorthogonally paired with respect to the Hilbert space inner product 



f (±,i) )Jf ( T>fe) \ =CjSjkf 

where C j = I ^(A^A'^A,) (W[cp+(x; A,), <p-(x; A,)]) 2 . 

Notice that the existence of a biorthogonal pairing automatically guarantees that the 
family of squared eigenfunctions evaluated at the double points of the spectrum of q form a 
linearly independent set. 

7. Appendix: Analytic Dependence of the Potential on the Deformation 

Parameter 

The reconstruction of a finite-gap potential of genus g relies on choice of 2g + 2 branch 
points in the complex plane and also on a choice of an effective divisor T> of degree g + 1 
(satisfying a certain reality condition) on the hyperelliptic Riemann surface defined by the 
branch points. Thus, to specify an isoperiodic deformation of the NLS potential q we must 
specify not only a deformation of the branch points but also a deformation of the divisor 
(or an equivalent set of auxiliary data). In this appendix, we show that, at the same time 
that a double point is being opened up by homotopic deformation, the auxiliary data given 
by the Dirichlet spectrum may be deformed in a way that guarantees that q is analytic in 
the deformation parameter e. This result is necessary for us to carry out the perturbation 
expansions in JUand 
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7.1. Finite-Gap Solutions and the Dirichlet Spectrum. The Dirichlet spectrum of an 
NLS solution q(x,t) which is L-periodic in x is defined as follows (see [TJ, [21])- Given a 
fundamental matrix solution $ for the AKNS system 

£)* ^K 1 < 33 > 

we construct the transfer matrix 

M(x, t; A) = $(x + L, t; \)$(x, t; A) -1 , 

There are two sets of Dirichlet eigenvalues: 

• The Dirichlet eigenvalues fik{x,t) are the A- values for which 

Mil - M 22 + Mi 2 - M 21 = 0. (34) 

Equivalently, these are the values for which there is a nontrivial solution of (|33|) 
satisfying the boundary condition (f>i + (f) 2 = at (x, t) and at (x + L,t). 

• The Dirichlet eigenvalues Vk{x,t) are the A- values for which 

M u - M 22 - i(M 12 + M 21 ) = 0. 

Equivalently, these are the values for which there is a nontrivial solution of (1331) 
satisfying i0 x + 2 = at (x, t) and (x + L,t). Also, these are the /i-eigenvalues that 
result when q is replaced by iq. 

The iik and v\. are dependent on x and t, satisfying a system of ODEs first derived by 
Ablowitz and Ma pp. However, for finite-gap solutions, all but finitely many of each set 
are locked to the double points of the Floquet spectrum. (This is true more generally for 
periodic NLS potentials whose Sobolev norm over one period is bounded; see |24j.) 

In order to get the equations defining the Dirichlet eigenvalues, we need to use the formulas 
for finite-gap NLS solutions and their associated Baker eigenf unction. We will briefly review 
their construction; further details are available in the book [2] or our paper [6]. 

One starts with a hyperelliptic Riemann surface S defined by 

2g+2 

i=i 

where the branch points Aj occur in complex conjugate pairs. This surface is compactified 
by the points oo±, where X g+1 /( — > ±1. After selecting a suitable homology basis cij,bj 
(see Figure [5]), one constructs holomorphic differentials Uj such that § ak ^j = 2iri5jk and 
meromorphic differentials df2i, df2 2 , df23 with prescribed singularities at oo ± and zero a- 
periods. Then 

/ \ „ , ^ , T ^(iVx + iWt-D-r) . . 

q ( x ,t) = A e M-^ + >Nt) L_____i , (36 ) 

where 

• 9 : C g — > C is the Riemann theta function, with period lattice generated by the 
vectors 2iriej (where {e^} is the standard basis of unit vectors) and with quasiperiods 
given by the columns of the Riemann matrix B defined by = <j> b uj; 

• V, W,and — r are the vectors of 6-periods of dOi, dfi 2 , df23 respectively; 



28 A. CALINI AND T. IVEY 

• the Abelian integrals n 1 ,0 2) ^3, with the branch point A 29+2 in the lower half-plane 
as basepoint, have asymptotic behavior 

ni(p)~±(A-| + o(i: 

exp(n 3 (P))~±r|A + (l) 

as P — > oo±, where A is real and positive. (Because they have nonzero 6-periods, the 
integrals are not well-defined on £, but are well-defined on the surface So obtained 
after cutting £ along the homology basis.) 

• The vector D is defined so that the zero divisor of 9(A(P) — D) is the positive divisor 
T> + linearly equivalent to T> — oo + , where A is the Abel map defined by 



A(P)= ! w 

J OO — 



and (V is the vector of holomorphic differentials aA (The reality condition on D is 

that D have zero real part.) 
A vector solution of fl33|) is provided at each P G X by the Baker eigenfunction -0, with 
components 

, ) = exp (i(fli(P)g + n 2 (P)t)) 9{A{P) + iVx + iWt - D)fl(D) 
' ~ exp(i(Ex-Nt)/2) 9(iVx + iWt - D)9(A(P) - D) 

and 

. exp (i(n 1 (P)x + Q 2 (P)t) + Q 3 (P)) 9(A(P) + iVx + iWt - D - r)fl(D) 
^ 2 ^' } ~ 1 exp(i(JVt - £a;)/2) 9(iYx + iWt - T>)6(A(P) - D) ' 

Although .<4.(P) and fli(P) are not well-defined on E because of nontrivial periods, we make 
the convention that the paths of integration for A(P) and flj(P) differ by a fixed path from 
oo_ to A 2 g+2 in the cut surface So, and this makes the Baker eigenfunction well-defined. 
(The rightmost theta-factors in the numerator and denominator, which were unfortunately 
omitted in [6], are necessary for this.) 

Let l : (A, () i — > (A, —() be the sheet interchange involution of the Riemann surface. When 
P is not a branch point, the Baker eigenfunctions i/)(P) and i^(lP) are a basis for solutions 
of Thus, to compute the initial values /ifc(0,0) and z/ fc (0,0) of the Dirichlet spectrum, 

we evaluate tp when t — and a; = or x = L. Using the facts that LV/(2ir) and LE/(2ir) 
are integer-valued, we get 

^i(0,0) = l, ^(0,0) = -i/(P), 

^(L,0) = exp (iL Ul x {P) - | M , ^(L,0) = -i/(P)^i(L,0), 

where 

/(P)=exp (Si 3(P))^ P) - D - r) 



0(A(P)-D) 

So, ipi+ip2 — is satisfied at both ends if and only either sin(Lf2i(P)) = or /(P) = — i. This 
reflects the fact that every double point of the Floquet spectrum is a Dirichlet eigenvalue 
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which is constant in x and t. On the other hand, the function f(P) is well-defined and 
meromorphic on E, with pole divisor T> + + oo + , so there are g + 1 Dirichlet eigenvalues 
/ifc(0, 0) that are not locked to the double points, and are located at the images of the points 
where f(P) = — i under the projection tt : £ — > C onto the complex A-plane. Similarly, the 
z/fc(0, 0) are either locked to double points or are images of the g + 1 points where f(P) = 1 
when x = t = 0. 

7.2. Deformations maintaining D = 0. 

Proposition 7.1. When D = 0, the points where f(P) = 1 (and hence the locations of 
Uk(0,0)) are A 2ff+2 and the branch points in the upper half plane (excluding A 29+2y ). -At the 
other branch points, f(P) = — 1. 

Before sketching the proof, we note that in order for the //-eigenvalues to deform continu- 
ously when further double points are opened up, the basepoint \2 g +2 must change continu- 
ously. Thus, we cannot take the convention (which we used in [6]) that it be located at the 
rightmost branch point in the lower half plane. Instead, we will keep it located at the branch 
point in the lower half plane that was originally part of the plane wave spectrum, located 
on the negative imaginary axis. Because this basepoint must be in the interior of the cut 
Riemann surface, we must rearrange our homology basis so that its cycles do not enclose this 
basepoint. (The extra cycle qq in the following diagram is not part of the homology basis.) 




Figure 5. Homology basis adapted to isoperiodic deformations, for genus 
g = 4. Branch cuts extend between complex conjugate branch points; solid 
paths lie on the upper sheet, dashed on the lower sheet; basepoint for Abelian 
differentials fi 2 , fi 3 is marked as A 29+2 . The extra cycles and I are used 
in the proof below. 



Proof. Because A 2ff+2 is the basepoint for Q 3 and .4.(A 29+2 ) = r/2, /(A 29+2 ) = 1. Let A 2g+ i = 
A 2s+2 in the upper half plane, and fix an integration path from A 2s+2 to A 2s+ i on the upper 
sheet on the right side of the branch cut. Because is a linear combination of differentials 
of the form (\ h /()d\, then f2 3 (A 29+ i) = — | § ao dfl 3 . Furthermore, because a + ai + . . . + a g 
is homologous to —I in the surface £ with the points oo± removed, then using the residue 
of df2 3 at oo + we get fis(A 2g+ i) = \ j^df2 3 = 7ri. Similarly, 

A(A 2g +i) = «4.(A 2g+2 ) - - / u> = + 7ril. 
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where uj is the vector of normalized holomorphic differentials u) k , and 1 is the vector whose 
entries are all equal to 1. Then, using the even-ness and periodicity of 9, 

f(\ \ p „ nr . ^(A(A 2g+1 ) - r) flfril - |r) g(|r-7ril) 

j(A 29 +i) = exp(vri)— — — — = - 1 = - — — — = -1. 

6'(w4(A2 S +i)) 0(7Til + ±r) 6»(7ril + ±r) 

Suppose Re(A 2 fe-i) < Re(A 2 c,+2) and lm(A 2 /t-i) > 0. Let b k be the cycle that runs from 
A 2 fc-i to A 2g +i along the upper sheet, and back along the same path along the lower sheet, 
whose projection to the A-plane encircles no other branch points. Then b k ~ bk + Yl a m, 
where the sum is over the a-cycles around branch cuts whose real part is strictly between 
that of A 2 fc_i and A 2g+ i. Letting the path of integration from A 2g +2 to A 2g +i be continued by 
the part of bk on the lower sheet, 

^3(A 2 fc-l) = f2 3 (A 2ff+ l) + | f dfi 3 = 7ri - \r k 

A(X 2k -i) = A(X 2g+1 ) + \l dVl 3 = §(r + B k ) + 7ri(l + ^e m ), 

where B k indicates the fcth column of the Riemann matrix. Then using evenness and peri- 
odicity of 9, 

6(A(\ 2k -i) - r) _ 0(f (-r + B k ) + 7ri(l + £ e m )) _ 9{\{v - B k ) + 7ri(l + £ e m )) 



and 



9(A(X 2k - 1 )) 0(i(r + J B fe )+7ri(l + Ee m )) 0(±(r - B k ) + 7ri(l + £ e m ) + B fc ) ' 
Using the quasiperiodicity property 6(z + 5^) = exp(— z k — ^B kk )9(z), 

f(x s = exp(7ri - |r fc )0(|(r - B k ) + 7ri(l + £ e m )) = 

1 2k ~ l} eM-Wk-B kk )+T:i-\B kk )9{\{v-B k )+m{l + ^ m )) 



Let A 2 fc = A 2 fc_i in the lower half-plane, and let the path from A 29+2 to A 2 fc_i be continued by 
the part of a k to the left of the branch cut. Then ^(A^) = fi 3 (A 2 fc_i) — | § df2 3 = fi 3 (A 2/ t_i) 
and A(A 2 fc) = A(\2k-i) — ^ e k- Using evenness and periodicity, 

9(A(X 2k ) - r) _ 0(§(-r + B k ) + 7ri(l - e k + £ e m )) _ 0(§(r - B k ) + 7ri(l - e fc + £ e m )) 



0(A(A 2fc )) " 0(i(r + Bft)+7ri(l-e fc + £em)) 0(|( r _ Sfc ) +7ri (i _ e ^ + £e m ) + 5* 
Then, using quasiperiodicity, 

^ . _ exp(?ri - |r fc )0(|(r - + 7ri(l - e fc + £ e m )) _ _^ 

exp(-|(r fe - B kk ) - \B kk )9(\{r - B k ) + 7ri(l - e k + J2 e m)) 



Similarly, assuming that Re(A 2 fc-i) > Re(A 2s+ i), lm(A 2 fe_i) > and \2 k = A 2 fe_i, we 
calculate that /(A 2 fc_i) = 1 and /(A 2 fc) = 1. □ 

Proposition 17.11 implies that, if we maintain the choice D = 0, any deformation of the 
branch points that is analytic in parameter e will induce an analytic deformation of the 
initial values u k (0, 0). This is also true if we increase the genus by analytically splitting a 
double point of the Floquet spectrum at e = into two branch points for e ^ 0, since the 
extra ^-eigenvalue that is unlocked at the higher genus still belongs to the Dirichlet spectrum 
as it limits to a double point when e = 0. It is also true for the isoperiodic deformations 
described in §2] that any other double point deforms analytically in e, and so any locked 
Dirichlet eigenvalues are analytic in e. 
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The /^-eigenvalues also deform analytically in e. For, when D = 0, the function / satisfies 
f(tP) = 1//(P). Thus, 

h(\) = f(P)+f(P)-\ A = vr(P) 

is a well-defined meromorphic function on the complex plane with pole divisor tc(D + + oo + ). 
The zeros of this function are fik(0, 0), so it has the form 

nLo(A-Mo,o)) 



h(X) 



r(A) 



where r(A) is a polynomial of degree g. At each branch point h(X) = ±2; the ^(0,0) are 
zeros of h(X) — 2, while the remaining branch points are zeros of h(X) + 2. Because the 
coefficients of the polynomials in the numerator of h(X) ± 2 must be analytic in e, it follows 
that the coefficients of the numerator of h(X) are analytic in e. Then the fik{0, 0) are analytic 
in e so long as they are distinct; however, distinctness holds for sufficiently small e at each 
deformation step. 

7.3. Trace Formulas. The NLS solution q(x,t) is determined by the Dirichlet spectrum 
and branch points by the trace formula^ 

9 

q(x, t) + r(x, t)=) j (X 2 k + X 2 k-i - 2fi k (x, t)) 

k=0 

9 

q(x, t) - r(x, t) = -i^J (^2k + X 2 k-i - 2u k (x, t)) , 

fc=0 

where we take the specialization r = — q for the focusing case. 

Given the branch points and the initial values fik(0, 0) and ^(0,0), the values of fik(x,t) 
and Uk(x,t) are determined by ODE systems derived by Ablowitz and Ma [TJ: 



dfi k _ ( C -E j¥fc ^) 



dx 
dt 



2fi, 



dfi k 2iC(/i fc ) 



+ 




+ r 



9 

3=0 



3=0 



where c = Em=i ^m, and 



duk 
dx 

di> k 
dt 



E 



2iCK) 



dvk 
ox 




+ 



5-E 



3=0 



3=0 



In these systems, the branch points X m and initial values /ifc(0, 0), z/ fc (0, 0) depend analyt- 
ically on e. Our choice D = implies that the initial values ^(0,0) are located at branch 
points, where ( = 0. Thus, the are constant in x and t, and are automatically analytic in 



4 These are adapted from [21], after making the changes q 1— ► — 2g, r 
to make their version of NLS and its Lax pair coincide with ours. 



-2r, A h- > — A which are necessary 
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In the system for the C(^fc) i s to evaluated along the upper sheet. (According to [TT] , 
the //-eigenvalues travel along the a-cycles of the Riemann surface.) Note that when e = 0, 
the initial values of one of the /i-eigenvalues is a double point, the limit of two complex 
conjugate branch points that coalesce as e — > 0. The above characterization of the /ij-(0, 0) 
as points where f(P) = — i, coupled with the fact that f(rP) = —f(P) ensures that the 
/ijt(0, 0) lie along the real axis. Thus, each £(/ifc(0, 0)) will be analytic in e, including at e = 0. 

It follows that the fj,k(x,t) and i/}.(x,t) are analytic in e for any x and t. Then the trace 
formulas imply that q(x, t) is analytic in e. 
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